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an  extensive  vehicle  sizing  capability  and  a state-of-the-art  optimization 
capability.  The  computer  memory  requirement  depends  upon  the  application 
of  the  system.  GTS  applications  on  oOOO  series  computers  typically  require 
125K  to  160K  octal  words  of  memory.  For  the  7UO0  series  computers,  GTS 
applications  typically  require  ldOK  to  L.OK  octal  words  of  small  core 
memory  ( S C M i and  50K  to  100K  octal  words  of  large  core  memorv  (LCMh 


This  volume  provides  the  potential  user  an  overview  of  the  G1S  system. 
Included  is  a summary  of  the  major  operational  capabilities  and  the 
structural  design  of  the  GTS  system.  V 
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PREFACE 


This  volume,  the  first  of  five  volumes  that  describe  the  Generalized 
Trajectory  Simulation  (GTS)  system,  provides  the  potential  user  with  an 
overview  of  GTS,  including  a summary  of  the  major  operational  capabilities 
and  structural  design  of  the  GTS  system.  The  remaining  volumes  are: 

i \ Volume  II:  GTS  Usage  Guide.  This  document  serves  as  a general 

^ ' usage  guide  for  GTS  and  includes  a set  of  example  problems,  a 

comprehensive  description  of  the  Generalized  Trajectory  Language, 
and  a discussion  of  the  trajectory  simulation  control.  In  addition, 
a set  of  appendices  contain  a master  reference  list  for  all  volumes 
and  supplementary  information  to  aid  tne  user  in  defining  his  problem. 

Volume  III:  GTS  Flight  Dynamic  Models.  This  publication  deals  v.  1 th 

the  GTS  library  of  flight  mechanics  and  flight  dynamics  models 
, utilized  for  trajectory  simulations. 

Volume  IV:  GTS  Numerical  Operators.  This  report  concerns 

J the  GTS  library  of  numerical  operators,  including  integration, 
optimization,  and  interpolation  operators. 

Volume  V : GTS  Weight  Estimation  Models  for  S .zing  Applications. 

This  volume  documents  the  GTS  library  of  weight  estimation 
models  utilized  'or  sizing  applications. 

This  report  was  prepared  by  D.  S.  Meder  and  J.  L.  Searcy. 
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.SECTION  1 
INI  RODUCTION 

Tills  document  is  the  first  of  a set  of  five  volumes  that  describe  the 
Generalized  Trajectory  Simulation  (GTS)  system.  The  intent  of  the  first 
volume  is  to  provide  the  potential  user  with  an  overview  of  the  capabilities 
and  structural  design  of  the  GTS  system.  Successive  volumes  provide  a 
more  detailed  description  of  the  various  components  of  the  GTS  system 
and  a description  of  the  input  data  required  to  execute  the  software. 

The  remaining  volumes  of  the  documentation  series  are: 

Volume  II.  GTS  Usage  Guide 

Volume  III.  GTS  Flight  Dynamics  Models 

Volume  IV.  GTS  Numerical  Operators 

Volume  V.  GTS  Weight  Estimation  Models  lor  Sizing  Applications 

The  Generalized  Trajectory  Simulation  (GTS)  system  comprises  a 
linked  set  of  digital  computer  programs  used  for  many  Aerospace 
Corporation  boost,  satellite,  and  reentry  vehicle /mi s sion  simulation 
studies.  GTS  satisfies  the  requirements  for  trajectory  simulations, 
vehicle  sizing,  and  optimization  capabilities  beyond  those  available 
in  othe  r computational  tools.  Architecturally,  the  software  system 
exploits  the  logical  concept  of  complete  separation  of  mathematical 
operations,  physical  modeling,  and  software  control  functions.  The 
system  software  is  completely  programmed,  with  some  minor  exceptions, 
in  a FORTRAN  language  compatible  with  CDC  6000/7000  series  computers. 
•Applications  may  be  submitted  to  the  system  via  the  normal  card  reader 
input  queue  or,  alternately,  from  remote  terminals. 

The  primary  design  consideration  for  the  GTS  system  has  been  the 
provision  of  an  easy-to-use  software  tool  that  permits  a variety  of  users, 
with  diverse  applications,  to  individually  balance  their  requirements  for 
completeness  and  accuracy'  with  cost  and  time  constraints.  Consequently, 
user -oriented  input  data  specification,  computational  efficiency,  wide 
applicability,  and  straightforward  software  development  and  modification 
features  have  been  primary'  objectives  ir.  the  design  of  the  program.  These 
objectives  have  been  accomplished,  in  part,  by  carefully  separating 


the  mathematical  operators,  the  engineering  models,  and  the  executive 
propram  control  functions.  The  logic  governing  the  interconnection 
between  program  components  is  under  input  processing  control.  As  a 
result  of  this  design,  the  user's  input  data  is  processed  and  related  ">  a 
directory  of  GTS  subroutines  so  that  the  minimal  software  configuration 
required  for  the  particular  application  is  automatically  loaded  and 
executed. 

The  modular  design  of  GTS  provides  an  exceptionally  convenient 
means  by  which  modifications  or  extensions  of  the  existing  capabilities 
can  be  implemented,  validation  and  verification  of  new  models  accomplished , 
and  overall  operating  efficiency  maintained  at  a high  level.  These  features 
are  of  particular  value  in  those  applications  which  are  exploratory  in  nature. 
The  major  components  of  the  GTS  system  are  given  in  Table  1. 

1.  1 GTS  Operational  Capabilities 

The  major  functions  and  capabilities  of  the  GTS  system  are  summarized 
in  this  section  to  provide  a convenient  introductory  reference.  Each  of  the 
significant  capabilities  is  briefly  described  in  subsequent  paragraphs  of  this 
section. 

1.1.1  Trajectory  Simulation 

11. e GTS  trajectory  simulation  capability  provides  a basic  framework 
to  define  and  control  a diversity  of  vehicle  and  mission  simulations,  The 
user  may'  define  new  physical  models  or  select  models  from  an  existing 
large  library  of  flight  dynamic  models.  A complete  description  of  the  existing 
models  is  provided  in  Volume  III.  The  GTS  trajectory  control  feature  allows 
a wide  degree  of  flexibility  for  the  specification  of  required  flight  profiles 
and  mission  objectives.  General  areas  of  application  of  the  trajectory 
simulation  capability  are  indicated  in  Table  d . 
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TO  SOLVE  THE  PARTICULAR  PROBLEM 
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1,  1.2  Optimization 


The  GTS  system  includes  a general  purpose  nonlinear  constrained 
optimization  capability  that  can  be  utilized  independently  or  in  conjunction 
with  the  trajectory  simulation  and  vehicle  sizing  capabilities.  The 
optimization  capability  is  designed  to  optimize  (i.  e.  . maximize  or 
minimize)  an  objective  function,  fix),  where  x is  a vector,  subject  to 
a set  of  nonlinear  equality  and  inequality  constraints.  The  operator 
is  mechanized  to  utilize  analytic  partial  derivatives,  if  available; 
otherwise  numerical  partial  derivatives  are  computed.  The  optimiza- 
tion features  of  the  GTS  system  are  described,  in  detail,  in  Volume  IV. 

1.  1.  ;>  Vehicle  Sizing  Analysis 

The  GTS  vehicle  sizing  models  provide  a vehicle  sizing  capability 
for  space  and  missile  vehicles  that  utilize  solid  propellant  rocket  motors. 
Liquid  propelled  vehicles  can  also  he  accommodated  within  the  sizing 
procedure.  They,  however,  require  the  development  of  specialized 
weight  estimation  models  to  represent  particular  characteristics  of  the 
desired  liquid  propellant  system. 

The  sizing  models  estimate  performance  sensitive  component 
weights  to  determine  a vehicle  configuration  consistent  with  realistic 
geometry  and  mission  constraints. 

Intended  applications  include  preliminary  design  studies,  booster 
subsystem  trade-off  studies,  and  growth  studies  of  existing  systems 
including  the  analysis  of  advanced  propellant  technology  or  new  launching 
concepts.  The  existing  library  of  weight  estimation  models  is 
completely  documented,  including  information  on  their  use  and  interface 
with  other  capabilities  of  the  GTS  system,  in  Volume  V. 

1.  1.4  Generalized  Trajectory  Language  (GIL) 

The  input  to  all  of  the  operational  capabilities  of  the  GTS  system  is 
via  a hiuhe  :•  level  symbolic  input  language  especially  designed  to  accom- 


modate  trajectory  simulations,  vehicle  sizing,  and  optimization  applications. 
GTL  not  only  controls  the  input  of  data  to  the  program,  it  also  is  used  to 
select  and  build  the  actual  software  configuration  that  is  loaded  into  the 
computer  and  executed. 

GTL  provides  a free  form,  user-oriented  input  method  that  is  straight- 
forward to  learn  and  use.  This  approach  encourages  concentration  on  the 
input  requirements  of  the  engineering  application  with  little,  or  no,  concern 
for  the  intricate  details  of  the  software.  The  input  language  also  allows  a 
user  to  include  FORTRAN  subroutines  as  a portion  of  the  input  data  deck. 
These  routines  are  processed  at  program  execution  time  and  utilized  in  a 
fashion  similar  to  any  other  model  within  the  GTS  system.  The  GTL 
language  is  described  in  Volume  II. 
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SECTION  2 i 

GTS  SYSTEM  CHAR  ACTERISTICS  ' 

A description  of  the  characteristics  of  the  major  GTS  system 
capabilities  is  presented  in  this  section.  A detailed  description  of  these 
capabilities,  including  a discussion  of  the  input  data  required  to  apply 
these  capabilities  to  specific  problems,  is  presented  in  the  remaining 
volumes  of  the  GTS  documentation  series. 

2.  1 Trajectory  Applications 

The  trajectory  simulation  capability  of  the  GTS  system  is  applicable 
to  a variety  of  vehicle  configurations,  mission  trajectory  profiles,  and 
mission  objectives,  as  indicated  by  Table  2.  Moreover,  the  modular 
structure  of  the  GTS  system  enables  the  user  to  satisfy  his  particular 
requirements  for  efficiency  and  completeness.  The  basis  of  the  trajectory 
simulation  capability  is  a set  of  flight  dynamic  models.  The  organization 
of  these  models  typifies  the  modular  design  of  GTS.  This  organization 
is  discussed  in  the  next  section.  An  overview  of  the  manner  by  which  these 
individual  models  are  incorporated  into  a trajectory  simulation  is  then 
discussed. 

2.  1.  1 Model  Types  and  Models 

The  flight  dynamics  models  are  organized  into  functional  units 
referred  to  as  model  types.  Associated  with  each  model  type  is  a set  of 
models.  Each  of  the  models  associated  with  a specific  model  type  simulates 
a particular  function  or  element  of  the  flight  dynamics.  For  example,  an 
atmosphere  model  simulates  the  atmosphere,  or  a weight  model  computes 
the  vehicle  weight.  A complete  list  of  the  trajectory  simulation  model  types 
and  their  corresponding  simulation  function  is  presented  in  Table  3. 

Clearly,  there  is  no  unique  way  to  simulate  these  functions.  Moreover, 
different  applications  require  varying  levels  of  complexity  of  the  functions 
being  simulated.  Consequently,  a variety  of  models  are  available  within 
each  model  type.  For  example,  there  are  currently  7 different  atmospheric 
models  available.  That  is,  each  model  within  a model  type  simulates  the 
same  engineering  function  but  in  a different  or  possibly  more  sophisticated 


Table  3.  Trajectory  Simulation  Model  Types  in  GTS 


MODEL  TYPE 
MNEMONIC 

FUNCTION 

1 N IT 

INITIALIZE  VEHICLE  POSITION  AND  VELOCITY 

GRAV 

SIMULATE  ACCELERATION  DUE  TO  THE  ATTRACTION 
Of  THE  CENTRAL  BODY 

ATM 

SIMULATE  ATMOSPHERE 

WIND 

SIMULATE  WINDS 

WEIGHT 

COMPUTE  VEHICLE  WEIGHTS 

GUID 

COMPUTE  GUI OANCE  COMMANDS.  OR  SIMULATE 
AIRBORNE  COMPUTER 

CONTROL 

SIMULATE  AUTOPILOT  OR  OPEN  LOOP  STEERING 

ATTITUDE 

COMPUTE  BODY  ATTITUDE 

AEROf 

COMPUTE  AERODYNAMIC  FORCES 

AEROM 

COMPUTE  AERODYNAMIC  MOMENTS 

AERO 

COMPUTE  AERODYNAMIC  FORCES  AND  MOMENTS 

ACTUATOR 

SIMULATE  ACTUATORS 

ENGINE 

SIMULATE  ENGINE 

PROPF 

COMPUTE  PROPULSIVE  FORCE  VECTOR 

PROPM 

COMPUTE  PROPULSIVE  MOMENTS 

TEOM 

INTRODUCE  THE  TRANSLATIONAL  EQUATIONS  Of  MOTION 

REOM 

INTRODUCE  THE  ROTATIONAL  EQUATIONS  Of  MOTION 

SENSOR 

SIMULATE  SENSORS;  e.g.,  ACCELEROMETERS,  GTROS, 
PLATFORMS,  IMUS 

DER 

INTEGRATE  AUXILIARY  EQUATIONS  Of  INTEREST 
(e.g.,  velocity  ImwsI 

NONDER 

COMPUTE  AUXILIARY  QUANTITIES  Of  INTEREST 
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manner.  That,  the  user  is  able  to  select  the  model  that  realistically  and 
efficiently  represents  his  problem. 

For  a particular  trajectory  simulation,  the  software  configuration  that 
is  actually  executed  includes  only  those  model  types  and  associated  models 
which  the  user  specified  for  that  simulation.  At  user-specified  events  new 
data  may  be  introduced  into  the  simulation  by  changing  the  data  associated  with 
a given  model,  by  changing  the  model  associated  with  a given  model  type,  or 
by  introducing  new  model  types  and  associated  models. 

2.  1.2  Trajectory  Profile  .Specification 

Within  the  GTS  system,  trajectory  simulations  are  performed  by 
describing  the  nesired  trajectory  profile  as  a series  of  phases  separated  by 
events.  A trajectory  is  defined  to  be  the  entire  vehicle  state-time 
history  of  interest.  A discrete  time  point  along  a trajectory  is  referred  to  as 
an  event.  The  user  may  specify  conditions  which  define  events.  These  condi- 
tions usually  refer  to  an  instantaneous  state  of  the  vehicle  as  computed  by  the 
program.  Examples  of  event  criteria  are: 

(i)  flight  path  angle  equal  to  -15  degrees 

(ii)  vehicle  altitude  increasing  and  equal  to  300000  feet 

(iii)  maximum  value  of  dynamic  pressure. 

Several  actions  may  result  from  the  detection  of  an  event.  For  example,  new 
data  may  be  introduced  to  the  simulation  which  will  alter  the  course  of  the 
trajectory,  or  information  may  be  output  which  will  leave  the  course  of  the 
trajectory  unaltered, 

A phase  is  defined  to  be  a continuous  portion  of  a trajectory  between 
two  successive  event  occurrences.  Thus,  a phase  is  always  initiated  and 
terminated  by  an  event.  Events  are  of  particular  interest  to  the  user,  since 
it  is  b\  the  response  to  us-,  r -s pecifie d events  that  new  data  is  introduced  to 
the  simulation  or  that  the  trajectory  is  terminated.  Precise  information 
concerning  the  types  of  events  and  the  specification  of  event  criteria  is 
presented  in  Volume  II. 
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2.  2 Sizing  Applications 

The  GTS  system  includes  a set  of  weight  estimation  models  which, 
combined  with  the  GTS  optimization  and  trajectory  simulation  features, 
form  a generalized  sizing  capability  for  vehicles  utilizing  solid  propellant 
rocket  motors.  Although  ultimate  component  design  is  not  precluded, 
the  sizing  capability  is  normally  for  preliminary  design  studies  to 
estimate  component  weights  which  will  determine  a propulsion  system 
configuration  consistent  with  realistic  vehicle  geometry,  performance, 
and  mission  constraints.  Parametric  weight  scaling  equations  are 
utilized  for  sizing  each  major  rocket  component.  The  sizing  capability 
does  not  require  the  specification  of  reference  designs  prior  to 
generating  results  and  provides  for  the  inclusion  of  technology  changes. 
Results  are  valid  for  propulsion  system  weights  between  3,000  and 
2, 000, 000  pounds. 

The  sizing  capability  utilizes  a model  type-model  structure  similar 
in  nature  to  the  trajectory  simulation  capability.  It  may  be  executed  as  a 
stand-alone  application  or  in  conjunction  with  the  trajectory  simulation. 
For  example,  the  sizing  capability  may  be  executed  to  evaluate  a vehicle 
geometry,  weight,  and  propulsion  quantities.  Subsequently,  at  the  user's 
option,  vehicle  parameters  computed  by  the  sizing  application  are 
available  as  input  data  for  a trajectory  simulation.  The  optimization 
process  is  invoked  to  satisfy  the  required  vehicle  geometry  constraint. 

Note  that  the  aforementioned  application  of  the  optimization 
capability  does  not  require  any  trajectory  simulations.  After  the 
optimized  vehicle  configuration  has  been  determined,  a trajectory 
simulation  utilizing  the  optimized  vehicle  may  be  requested. 

A variation  of  the  procedure  just  described  is  the  application  of 
the  optimization  operator  to  the  satisfaction  of  coupled  vehicle /trajectory 
constraints.  In  this  instance,  the  objective  might  be  the  determination  of 
a set  of  vehicle  parameters  that  satisfy  a set  of  vehicle  geometry 
constraints  and  maximizes  the  range  of  the  vehicle.  For  such  an 
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application,  the  sizing  capability,  the  trajectory  simulation  capability 
and  the  optimization  capability  are  all  executed  as  a part  of  the  same 
application.  Additional  information  on  the  individual  sizing  models 
and  examples  of  the  use  of  the  sizing  capability  are  presented  in 
Volume  V. 

2.  3 Numerical  Operators 

The  extensive  GTS  library  of  physical  and  engineering  models  is 
supplemented  with  a compatible  set  of  numerical  operators.  These 
operators  are  d is tingui shea  from  tne  physical  and  engineering  models 
in  that  they  provide  methods  for  the  solution  of  mathematical  problems 
associated  with  scientific  and  engineering  simulations.  The  three 
major  categories  of  numerical  operators  available  within  the  GTS 
system  are;  (1)  integration  operators,  (2)  optimization  and  search 
operators,  and  (3)  interpolation  operators.  Analogous  to  the  engineering 
operators,  GTS  provides  a choice  of  operators  in  each  category. 
Furthermore,  the  numerical  operators  are  independent  of  the 
engineering  models  ir.  the  sense  that  any  of  the  numerical  operators 
may  be  specified  in  conjunction  with  any  particular  model  configuration 
the  user  has  selected.  For  example,  any  one  of  the  available  integration 
operators  nay  be  selected  for  a particular  trajectory  simulation. 
Similarly,  any  one  of  the  available  optimization  algorithms  may  be 
executed  as  a part  of  the  trajectory  optimization  capability  or  vehicle 
sizing  capability.  Consequently,  the  engineering  analyst  may  choose 
the  numerical  algorithms  which  provide  the  desired  balance  between 
accuracy  and  efficiency  for  his  particular  problem.  These  numerical 
algorithms,  along  with  the  associated  input  data,  are  described  in 
Volume  IV. 

2.  3.  1 Integration  Techniques 

While  the  integration  techniques  in  GTS  are  general  purpose,  these 
techniques  have  been  chosen  with  particular  attention  to  their  viability  in 
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solving  sets  of  differential  equations  associated  with  trajectory  simulations. 
Variations  of  mission  profiles,  accuracy  requirements,  and  computational 
efficiency  dictate  that  several  techniques  be  available.  These  techniques 
are  Runge-Kutta  fixed  step  (fourth  order),  Adams  - Moul  ton  fixed  step 
(rn-th  order,  1 < m £ 8).  and  A da  ms-Moul  ton  variable  step,  (m-th  order, 

1 £ m £8).  The  integration  operators  are  not  restricted  to  a pre-set 
number  of  differential  equations.  The  number  of  differential  equations  is 
dependent  upon  the  particular  simulation,  and  typically  varies  from  as  few 
as  six  to  over  one  hundred  equations. 

2.3.2  Optimization  Capability 

The  GTS  system  contains  two  state-of-the-art  nonlinear  constrained 
parameter  optimization  operators.  They  are  designed  to  solve  the  following 
problem.  Determine  the  values  of  the  n parameters,  or  independent 

variables,  x^,  such  that  the  scalar  function  f(x)  = f(Xj,  x^»  . ••  , x^) 

is  optimized  (maximized  or  minimized),  subject  to  equality  constraints  of 
the  form 

c.  (x)  = 0 , i = 1 , 2,  . . . , k 

and  inequality  constraints  of  the  form 

c.(x)  >0  , i = k+1 , . . . , m 

Special  cases  of  the  above  problem  include  constraint  solving  problems  in 
which  there  is  no  objective  function  and  the  number  of  variaules  is  greater 
than  or  equal  to  the  number  of  constraints,  and  nonlinear  least  squares 
problems  in  which  there  is  no  objective  function  and  the  number  of 
constraints  is  greater  than  the  number  of  variables.  Also,  optimal  control 
problems  may  be  solved  using  the  optimization  techniques  available  in  GTS. 

These  optimization  and  search  techniques  are  fully  incorporated  into 
the  GTS  system.  The  complete  library  of  physical  models  is  available  as 
part  of  the  definition  of  an  optimization  problem.  The  input  required  to 
define  optimization  problems  is  compatible  with  the  input  required  to  define 
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trajectory  simulation  and  vehicle  sizing  problems.  Thus,  the  definition 
of  the  optimization  problem  does  not  require  a separate  specification  of 
the  associated  trajectory  simulation  or  weight  estimation  problem.  To 
illustrate,  assume  that  a user  has  defined  a trajectory  simulation,  and  as 
a part  of  further  analysis,  he  wishes  to  define  an  associated  trajectory 
optimization  problem.  The  only  additional  input  required  is  a specification 
of  which  trajectory  input  quantities  are  to  be  the  independent  variables  and 
which  computed  quantities  are  to  be  the  objective  function  and/or  constraints. 

In  this  regard,  any  input  quantity  may  be  an  optimization  variable.  Thus, 
for  example,  an  independent  variable  may  be  a tabular  value,  a model  input 
value,  or  a quantity  which  is  part  of  an  event  specification.  Likewise, 
the  user  has  complete  freedom  in  choosing  which  quantities  form  the 
objective  function  and  constraints. 

Furthermore,  the  independence  o'  the  physical  models  and  optimization 
provides  the  user  the  capability  to  formulate  many  diverse  types  of  problems. 
As  indicated,  a variety  of  trajectory  optimization  problems  may  be  formulated, 
such  as  boost  trajectory  problems,  reentry  problems,  or  orbital  transfer 
problems.  Likewise,  vehicle  sizing  applications  may  be  defined  which  do 
not  require  trajectory  simulations,  or  vehicle  sizing  applications  may  be 
specified  which  require  both  the  vehicle  sizing  and  trajectory  simulation 
capabilities.  Finally,  the  operators  may  also  be  used  to  obtain  the  optimal 
solution  to  a large  class  of  problems  having  similar  mathematical  character- 
istics to  trajectory  and  vehicle  sizing  problems. 

Some  of  the  optimization  algorithms  utilize  partial  derivative  informa- 
tion. If  the  information  can  be  computed  analytically,  the  optimization 
operator  will  effectively  use  this  information.  Otherwise,  a scheme  for 
obtaining  numerical  partials  is  implemented. 

2.  4 General  Applicability  of  the  GTS  System 

Many  of  the  capabilities  of  the  GTS  system  have  applicability  to  areas 
other  than  trajectory  simulation  or  vehicle  design.  The  numerical  operators, 
for  instance,  provide  numerical  methods  for  solving  mathematical  problems 
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which  may  arise  independent  of  any  vehicle  design  or  trajectory  simulation 
considerations.  For  example,  the  integration  operators  provide  a capability 
for  the  solution  of  any  set  of  first  order,  initial  value  differential  equations. 
Likewise,  the  optimization  operators  are  designed  to  solve  a general  non- 
linear parameter  optimization  problem.  The  flexibility  of  the  GTS  input 
characteristics  permits  these  operators  to  be  properly  invoked  to  solve 
the  particular  problem  of  interest.  This  capability  is  illustrated  in 
Volume  II. 

2,5  Generalized  Trajectory  Language  (GTL) 

GTL  is  a higher  level  input  language  included  as  an  integral  part  of 
the  GTS  system.  GTL  interprets  card  images  and  formulates  all  the 
input  into  a form  acceptable  to  the  internal  GTS  input  processor.  GTL 
provides  a use r -oriented  format  that  permits  a user  to  define  the  application 
(trajectory,  vehicle  sizing,  optimization,  etc.), the  individual  models 
required,  the  associated  numerical  operators,  and  all  of  the  required  input 
data  values  in  a free  form  mnemonic  oriented  fashion.  A complete  descrip- 
tion of  GTL  is  part  of  Volume  II. 


To  illustrate  GTL,  consider  the  following  examples  of  GTL  data 
specification.  Relative  to  the  trajectory  simulation  capability,  the  following 
examples  illustrate  the  specifications  of  events  in  GTL.  Assume  an  event 
occurs  when  the  criterion  that  the  variable  VR  is  equal  to  1000  is  satisfied. 
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The  GTL  s pec  if  ication , that  is  the  actual  input  that  users  must  supply,  for 
this  event  is 


DBQOE1  SDQEin!]  HBlBIEffinniniDE!]  3l3BH03B]05]E!]nJE0EBJ00E 

■■  Kiel 


Next,  assume  an  event  occurs  when  the  criterion  that  the  variable  HFT 
is  equal  to  the  variable  HATM,  and  HFT  is  increasing  (point  A on  the  graph). 


HATM 


1 he  GTL  specification  for  this  event  is 


I naans  ehwoiei DCBEiEoiBnurnEE]  jp^SBigjig 

S 


31  32  33  34 


mm 


FT[  + 


Alternatively,  if  an  event  is  to  occur  when  HFT  is  equal  to  HATM  and  HFT 
is  decreasing  (point  B on  the  graph),  then  the  GTL  specification  of  this  event  is 


CSUBN 


A common  method  for  specifying  functional  relationships  in  GTS  is 
via  tabular  data.  GTL  provides  a convenient  format  for  tabular  input.  As  an 
illustration,  assume  that  the  quantity  CSUBN  is  a function  of  two  independent 
variables,  ALPHAT  and  MACHNO,  as  represented  in  the  following  graph. 


The  cor  re sponding  GTL  tabular  input,  assuming  quadratic  interpolation, 
is  the  following. 


Iii  order  to  provide  a convenient  method  for  executing  the  optimization 
capability,  GTL  includes  an  optimization  input  format.  The  format  provides  the 
user  with  an  easy-to-use  mechanism  for  specifying  optimization  input  data. 

For  example,  the  choice  of  the  optimization  operator  is  specified  by  a 
statement  of  the  form 


1L22  23124  25 C 6 2 7(28  29  30  31,32,33  34  35  34 '37- 
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A constraint  that  the  variable  VR,  evaluated  at  event  PO70,  be  greater 
than  or  equal  to  1000  is  specified  by 


1*  19  20 


(GDlV  GLUPD.  f.  30. 


An  additional  GTL  capability  allows  user  coded  FORTRAN  routines 
to  be  included  in  the  GTL  input  data.  These  routines  are  automatically 
inserted  into  the  GTS  simulation.  Thus,  the  user  can  request  that 
unique  auxiliary  computations  be  made  or  that  minor  modifications  to 
existing  GTS  models  be  made  without  modifying  the  existing  GTS  models. 


z.  o Input/Output  Gapaoiiitios 

The  GTS  system  allows  the  user  to  specify  his  input  data  and  receive 
his  output  data  in  v irious  forms  and  formats.  The  primary'  means  of 
specifying  input  data  is  by  means  of  punched  cards.  However,  the  user  may 
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access  input  data  which  resides  on  a disk  file  or  magnetic  tape.  Special 
purpose  input  data,  such  as  radar  tracking  data,  may  be  input  to  a problem. 
The  system  requires  the  user  to  describe  the  format  ol  the  input  data.  That 
is,  the  number  ol  words  per  data  frame  or  the  relative  locations  of  the 
particular  data  items  of  interest  are  specified  to  the  system  via  GTL  input. 

GTS  has  the  capability  of  producing  several  types  of  output  all  under 
the  control  of  the  user.  Included  are  printed  output  in  a standard  format, 
printed  output  and/or  case  summary  print  which  is  formatted  by  the  user, 
binary  data  files  for  subsequent  use  in  other  software  tools,  and  graphical 
output  processors. 

The  standard  output  print  (block  print)  follows  the  model  types  of  GTL 
input  data.  That  is,  the  printed  data  is  formed  into  distinct  blocks  where 
each  block  is  associated  with  a particular  model  type  (e.g.  , AEROF, 
CONTROL).  A data  item  may  be  a floating  point  number,  an  integer 
number,  a logical  value,  or  an  alpha-numeric  character  string.  The  user 
may  request  that  blocks  be  edited  to  reduce  the  output  or  that  data  items 
which  are  not  included  in  a standard  print  block  be  printed  in  a supplementary 
data  block.  A block  print  key  is  printed  at  the  beginning  of  each  case. 

User-defined  output  is  handled  by  a special  output  operator.  Any  GTS 
variable  or  user-defined  variable  name  is  an  acceptable  output  quantity. 
Multiple  executions  of  the  output  operator  are  permissible,  and  each  call 
may  reference  a different  format  statement  and  list  of  variables.  In  addition, 
if  logical  tests  are  required  to  determine  what  to  print  or  when  to  print,  a 
subroutine  may  be  written  and  input  with  the  GTL  input  data  to  perform  these 
functions. 

Plots  are  obtained  by  writing  the  output  data  to  be  plotted  on  files 
using  the  special  output  model.  The  plotting  capabilities  are  executed  as 
pos t- pr oce s s o r s and  are  accessed  by  means  of  GTL  input.  Also  available 
is  a world  map  plotting  capability.  The  user  may  also  request  that  any  or 
all  input  data  tables  (interpolation  or  step  tables)  be  plotted. 
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2.  7 Trajectory  Simulation  Example 

As  an  illustration  of  the  use  of  the  GTS  system,  consider  the  simulation 
of  the  three  deg  ree -of -freedom  trajectory  for  a satellite  in  orbit.  The  data 
specified  to  define  the  orbit  and  establish  the  initial  position  of  the  vehicle  are: 


apogee  altitude 

- 

3000  nautical  miles 

perigee  altitude 

= 

1000  nautical  miles 

inclination  of  orbit 

= 

20° 

argument  of  perigee 

- 

85° 

t rue  anomaly 

- 

ro 

a- 

o 

o 

satellite  longitude 

- 

-119° 

A spherical  non-rotating  earth  model  is  used  in  this  example.  The 
integration  method  is  8th  order,  variable  step,  Adams -Moulton,  The  initial 
integration  interval  is  8 seconds.  Trajectory  data  is  to  be  printed  at  apogee 
and  perigee  (i.  e.  , whenever  the  inertial  flight  Dath  angle  is  zero). 

The  initial  configuration  of  the  orbit  is  shown  in  Figure  1.  The 
data  as  it  appears  on  a GTS  coding  form  is  shown  in  Figure  2.  A computer 
generated  listing  of  the  CTL  input  and  trajectory  output  is  also  included. 
Additional  examples  which  illustrate  the  capabilities  of  the  GTS  system  are 
detailed  in  Volume  II. 
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The  preceding  BLOCK  PRINT  KEY  displays  print  key  lines  which 
contain  mnemonic  program  symbols.  The  print  key  lines  are  identified  by 
line  numbers.  These  line  numbers  are  used  to  relate  the  lines  of  printed 
output  in  the  BLOCK  PRINT  with  the  print  key  lines  contained  within  the 
BLOCK  PRINT  KEY. 

Brief  definitions  of  the  output  quantities  are  presented  below. 

The  output  quantities  are  identified  by  their  mnemonic  program  symbol,  and 
are  ordered  by  line  number. 


LINE  1 

TIME 

HFT 

GCLV 

GDLV 

LV 

W 


Current  time  within  the  dynamic  system 
Altitude  of  the  vehicle  (feet) 

Geocentric  latitude  of  the  vehicle 
Geodetic  latitude  of  the  vehicle 
Longitude  of  the  vehicle 
Total  vehicle  weight 


LINE  2 

DTNOM 

H 

RV 

ETATOT 


Nominal  integration  step  size 
Altitude  of  the  vehicle  (nautical  miles) 

Radial  distance  from  center  of  earth  to  vehicle 
Total  load  factor 


DA- 

LINE  3 
VI 

GAMMA! 

AZI 

VR 

GAMMAR 
A ZR 

LINE  4 

XECI 

YECI 

ZECI 

DXECJ 

DYECI 

DZECI 


Total  weight  flow-rate  from  vehicle 


Magnitude  of  the  inertial  velocity  vector 
Inertial  flight  path  angle 
Inertial  azimuth 

Magnitude  of  the  earth  relative  velocity  vector 
Earth  relative  flight  path  angle 
Earth  relative  azimuth 


The  ECI  X coordinate  of  vehicle 
The  ECI  Y coordinate  of  vehicle 
The  ECI  Z coordinate  of  vehicle 
Time  derivative  of  XECI 
Time  derivative  of  YECI 
Time  derivative  of  ZECI 


pos  ition 
position 
pos  ition 
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UNE  5 


DDXECI 

DDYECI 

DDZECI 

SDDXECI 

SDDYECI 

SDDZECI 


Time  derivative  of  DXECI 
Time  derivative  of  DYECI 
Time  derivative  of  DZECI 
The  ECI  X component  of  sensed 
The  ECI  Y component  of  sensed 
The  ECI  Z component  of  sensed 


accele  ration 
accele  ration 
accele  ration 


LINE  7 

GECI(l) 

GECI(2) 

GECI(3) 

ILV 


The  ECI  X component  of  gravitational  acceleration 
The  ECI  Y component  of  gravitational  acceleration 
The  ECI  Z component  of  gravitational  acceleration 
Inert;al  longitude  of  the  vehicle 


LINE  1 1 

RANGRLS 
RANGRIP 
R ANGILS 
RANGIIP 


Relative  range  from  launch  site  to  current  position 
Relative  range  from  initial  position  to  current  position 
Inertial  range  from  launch  site  to  current  position 
Inertial  range  from  initial  position  to  current  position 


i 


LINE  30 
ORBIT 

TYPE  Type  of  orbit 

SOVRCEA  Ratio  of  s emilatus  - rectum  to  equatorial  radius 
KE/PL  Ratio  of  kinetic  energy  to  potential  energy 


LINE  31 

INCNODE 

ECC 

DOMEGAP 

PERIOD 

DLNODE 

SMAORB 


Orbit  inclination  at  next  ascending  node 
Orbit  eccentricity 
Inertial  motion  of  apsides 
Orbit  pe  r iod 

Inertial  motion  of  ascending  node 
Orbit  3emimajor  axis 


LINE  32 


INC  LV 

TRANOMV 

OMEGA V 

TNODE 

LNODE 

MUNODE 


LINE  33 
HA 

VAPOG 

OMEGAA 

GDLVA 

LVA 

RA 


LINE  34 
HP 

VPERG 
OMEGA P 
GDLVP 
LV  P 
RP 


LINE  3 5 


TAPOG 

MUAPOG 

TPERG 

MUPERG 


Local  orbit  inclination 
T rue  anomaly 

Argument  of  latitude  of  current  position 
Time  of  next  ascending  node 
Relative  longitude  of  next  ascending  node 
Inertial  range  angle  from  present  position  to  next 
ascending  node 


Altitude  of  next  apogee 
Inertial  velocity  at  next  apogee 
Argument  of  latitude  of  next  apogee 
Geodetic  latitude  at  next  apogee 
Relative  longitude  at  next  apogee 

Radial  distance  from  the  center  of  the  earth  to  the 
next  apogee 


Altitude  at  next  perigee 

Inertial  velocity  at  next  perigee 

Argument  of  latitude  of  next  perigee 

Geodetic  latitude  at  next  perigee 

Relative  longitude  at  next  perigee 

Radial  distance  from  the  center  of  the  earth  to  the 

next  perigee 


Time  of  next  apogee 

Inertial  range  angle  from  present  position  to 

next  apogee 

Time  of  next  perigee 

Inertial  range  angle  from  present  position  to 
next  perigee 
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2.8  Remote  Terminal  Execution  of  GTS 

All  of  the  capabilities  of  the  GTS  system  are  available  for  remote 
terminal  use  in  conjunction  with  the  CDC  supplied  INTERCOM  feature. 

The  actual  methods  of  operation  will  depend  on  the  hardware  characteristics 
of  the  terminal  to  be  used. 

The  output  from  the  job  may  be  obtained  by  several  different  methods. 
The  user  may  request  the  print  file  to  be  directed  to  his  terminal  where  he 
may  either  list  the  entire  file  or  make  use  of  INTERCOM  editing  procedures 
to  print  items  of  particular  interest.  If  the  user  has  a line  printer  available, 
he  may  list  the  output  by  that  means.  Alternately,  the  output  can  be  directed 
to  the  central  site  printer. 

Figure  3 shows  a portion  of  the  input  required  to  execute  the  GTS 
system  from  a remote  terminal  via  the  INTERCOM  system.  The  INTERCOM 
system  directs  the  flow  of  information  and  data  between  a central  computer 
and  a number  of  remote  terminals  and  is  compatible  with  CDC  6000/7000 
computers.  For  this  example  the  trajectory  simulation  that  was  defined  in 
section  2.  7 is  executed. 

The  reader  can  easily  recognize  the  correlation  of  the  data  given  in 
Figure  3 and  the  data  given  in  Figure  2.  The  additional  symbols  and 
commands  are  specific  to  the  INTERCOM  system  and  are  explained  in  the 
INTERCOM  User's  Guide  available  from  Control  Data  Corporation. 
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The  followinq  statements  are  remote  terminal  inputs  to  GTS 


5l0=JUb  ORBIT  LCCLivlRIC  ORBIT 
520=bLCIN  EVENT  bt^ULNCL 
5 3G  = PO10 


54  0 = 

PhlO  INITIALIZE 

VEHICLE 

ON  ORBIT 

550  = 

IN  IT 

INI TM  30 

BAD 

= 3000 

560  = 

H PD 

= 1000 

570  = 

INCLNIN 

= 20 

5oG  = 

UP 

= £ 5 

590  = 

TANOH IN 

= 260 

600- 

LTFOl 

= -119 

610  = 

^RAV 

GRAVMl 

CBCON 

= 2 

620  = 
630  = 

TEGM 

TEOMMl 

FOTATL 

= 0 

640  = 

INTGRA 

INTGRMl 

METHOD 

= : A Ni  VARIABLE  STEP 

650  = 

ORDER 

= 8 

660  = 
670  = 

OUT 

OUTPUT 

DT  IN 

= 8 

6b0  = 9 

6 9 0 =R  E 2 0 GAMMA  I = 0 
700  = 9 

7 1 0 = POl 0 0 TIME  = 12000  /STOP/ 
7 2 0 = LND  EVENT  SEQUENCE 

7 30  = ? 

7 4 0 = ? 

7 5 0 = END  OF  CA.sE 


The  followinq  statements  are  INTERCOM  corrands  used  to  send 
the  above  inputs  to  the  central  computer. 


save ,orbi t , n 
catch, or  bit, input 


Fig.  3.  Remote  Terminal  Input 
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SECTION  3 

GTS  SYSTEM  ARCHITECTURE 

The  GTS  system  consists  of  a compatible  set  of  processors,  operators, 
and  models  that  can  be  t'Lexibly  selected  to  prepare  a minimuni  software 
configuration  customized  to  the  requirements  of  a particular  application. 

The  primary  functional  capabilities  (trajectory,  vehicle  sizing,  optimization) 
can  be  invoked  in  any  logical  hierarchal  fashion.  Each  of  these  capabilities 
invokes  and  loads  only  the  models  and  operators  required  to  accomplish  its 
purpose,  All  of  the  models  and  operators  are  interfaced  through  a 
sophisticated  data  structure  and  corresponding  set  of  data  operators.  This 
feature  permits  independently  developed  models  to  be  properly  interfaced, 
via  the  data  structure,  into  a coherent  simulation.  Consequently,  the 
system  support  structure  also  facilitates  the  extension  or  altering  of 
simulations  by  adding,  deleting,  or  changing  model  subroutines.  In 
conjunction  with  GTL,  the  system  will  load  and  execute,  in  a stand-alone 
fashion,  any  individual  engineering  model  subroutine  using  input  defined 
parameters  and  providing  corresponding  outputs.  This  feature  provides 
a powerful  tool  for  the  modular  validation  and  verification  of  large-scale 
simulations. 

The  complete  system  has  been  automated  to  eliminate  bothersome 
and  time-consuming  operational  overhead.  The  required  processors 
are  invoked,  in  the  proper  sequence,  on  the  command  of  C DC  control 
cards  that  are  automatically  generated  for  the  user  by  a control  card 
generator.  This  allows  the  user  to  concentrate  on  the  application  and 
not  on  the  physical  operation  of  the  software. 

The  component  members  of  the  GTS  system  architecture  are 
indicated  in  the  following  paragraphs  of  this  section  along  with  a brief 
description  of  their  function. 

3.  1 Control  Card  Generator 

The  control  card  generator  constructs  proper  sets  of  CDC  control  cards 
to  execute  the  GTS  system  for  a variety  of  applications.  Options  are 
available  to  invoke  features  of  a maintenance  system,  provide  for 
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modifications  and  additions  to  the  library  of  operators  and  models  and 
to  select,  load,  and  execute  a simulation. 

3.  2 GTL  Processor 

The  GTL  processor  scans  the  GTL.  input  data  to  construct  files  of 
parameter  and  variable  data  required  by  a given  application  plus  determining 
the  functional  processors,  associated  models,  and  program  structure 
required  for  the  application.  The  program  structure  information  is  passed 
to  a subsequent  pre -processor  that  selects  the  required  program  elements 
from  a library  and  constructs  a tailorized  version  of  GTS.  The  parameter 
and  variable  information  is  saved  on  temporary  files  for  use  in  the 
application  as  it  is  required. 

3.3  GTS  Pre  -proces  sor 

This  program  module  accepts  files  of  program  structure  information 
compiled  by  the  GTL  processor.  This  pre -proces sor  utilizes  this 
information  to  select  the  required  program  models  from  the  GTS  system 
library.  These  models  and  all  implicitly  required  supporting  subroutines 
are  structured  into  a computer  program  specific  to  the  given  application. 

The  computer  program  is  further  processed  by  the  CDC  system  loader 
and  made  available  for  execution.  This  procedure  eliminates  many  of  the 
problems  normally  associated  with  the  ever -expanding  computer  memory 
requirements.  The  size  and  number  of  operators  and  subroutine  models 
contained  in  the  GTS  library  is  allowed  to  grow'  in  an  open-ended  fashion. 
This  grow'th  does  not  result  in  an  ever  increasing  requirement  for  computer 
memory. 

3.  4 GTS  Executives  and  Model  Library 

The  primary  executive  operators  of  the  GTS  system  were  identified 
in  Section  1,  1.  These  functional  operators  (trajectory,  vehicle  sizing, 
optimization)  are  supported  by  an  extensive  library  of  numerical  operators 
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(integration,  interpolation,  etc.),  data  operators,  and  physical  and 
engineering  model  subroutines.  Libraries  are  maintained  in  FORTRAN 
source  language  and  relocatable  binary  format  for  CDC  o000/7300  series 
computers . 

3.  5 GTS  Post-Processors 

A series  of  post-processing  elements  are  available  to  the  user  of 
the  GTS  system.  Communication  with  these  processors  is  via  the  use  of 
disk  or  tape  files  generated  during  the  execution  of  a GTS  simulation. 

The  primary  post-processors  provide  a variety  of  graphical  output 
capabilities.  These  graphical  features  provide  the  ability  to  plot  any 
GTS  variable  versus  any  other  %'ariable.  Features  are  also  available  to 
display  information  on  p re-printed  or  computer  generated  world  map 
backg  rounds . 
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